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As a result of the development of wind farms, the gas — steam blocks, which shall quickly ensure energy supply in 


case the wind velocity is too low, are introduced to the energy system. To shorten the start-up time of the gas — 


steam and conventional blocks, the structure of the basic components of the blocks are changed, e.g. by reducing 


the diameter of the boiler, the thickness of its wall is also reduced. The attempts were also made to revise the cur- 


rently binding TRD 301 regulations, replacing them by the EN 12952-3 European Standard, to reduce the allowa- 


ble heating and cooling rates of thick walled boiler components. The basic assumption, on which the boiler regu- 


lations allowing to calculate the allowable temperature change rates of pressure components were based, was the 


quasi — steady state of the temperature field in the simple shaped component, such as a slab, cylindrical or spher- 


ical wall. 
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Introduction 


The quasi — steady state occurs in structural compo- 
nents of boilers and turbines during a long heating or 
cooling process running at a constant rate of change of 
temperature. In this paper, a more detailed analysis of the 
quasi — steady state was conducted for non — weakened 
walls and for walls weakened by openings. The use of 
formulas for circumferential stress at the edges of the 
openings to determine the allowable heating and cooling 
rates for the assumed pressure was also presented. Those 
formulas can be used in real life, as the analysis of the 
temperature and stress fields in the quasi — steady state 
can be quickly conducted with the use of the Finite 
Elements Method (FEM), even for components of a 
complex structure. 

The quasi — steady state can be characterized by the 


Nomenclature 
c specific heat of the metal, J/(kg:K) t 
di inner diameter of the drum, m T 


following factors: 
e constant temperature difference between the se- 
lected points of the component (Fig. 1) 

e constant, time independent, thermal stress (Fig. 2). 

The results presented in Figs. 1 and 2 were obtained 
for the following data: r;, = 0.85 m, r, = 0.94 m, c = 511 
J(kg:K), p = 7775 kg/m’, k = 47.3 W/(m-K), E = 201105 
MPa, £ = 1.21x10° 1/K, v= 0.288. 

Let's consider the simplifying assumption that the 
thermo-physical properties of the wall material: k ,c, p, E, 
and v are constant and determined for the average tem- 
perature over time and the wall thickness. The transient 
heat transfer equation in the quasi — steady state can be 
then simplified to form [1, 2] 


VT = (1) 
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time, s 
temperature, °C or K 


dwo inner diameter of the downcomer, m 

E modulus of elasticity, MPa On 

h heat transfer coefficient, W/(m*-K) ar 

k thermal conductivity, W/(m:K) B 
absolute pressure in the drum, MPa K 

Po atmospheric pressure, MPa v 


gauge pressure in the drum p, = p — po, 


Greek Symbols 


stress concentration factor for circumferential 
stress caused by pressure 

stress concentration factor for circumferential 
thermal stress 

linear thermal expansion coefficient, 1/K 
thermal diffusivity, m?/s 
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Fig. 1 The temperature difference between the internal and Fig.2 Circumferential stress on the inner and outer surface of 


the external surface of the cylinder during heating at a 
constant rate vr = 3 K/min. 


oT : 
where vr ap ee is the rate of wall tempera- 
t 


ture change. The symbol V? denotes the Laplace operator. 
The thickness of the wall is constant, and the external 
surface r = r, of the pressure component is thermally 
insulated 

dT 


dr r=, 


=0 (2) 


The temperature on the internal surface r = r;, rises 
linearly in time at a rate of vr 


T| Tn =Vrt (3) 


Expressing the Laplace operator in a cylindrical or 


Tri 


the boiler drum during heating at a constant rate: vr = 
3 K/min. 


spherical coordinate system leads to the specific form of 


the Eq. (1) 
1 (r2) (4) 
r” dr dr K 
where m = 1 for the cylindrical component and m = 2 
for the spherical component. After double integration of 
the above equation and after determining the constants 
from the boundary conditions (2) — (3) for the cylindrical 
wall, the following expression is obtained 


T =T,, 4 fl G re 2r? In 7 (5) 
4K r, 


The average temperature over the wall thickness is 


calculated from the formula 


m+ Rom 
En = el ot [ "Tar (6) 
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Substituting Eq. (5) into Eq. (6) with m = 1 and inte- 
grating gives the average temperature T, on the thickness 
of the cylindrical wall 

2 2 4 
T- patil -1)(3u aes Inu 
K 8(u? -1)(w-1) 
where: u =1,/r;, - the ratio of the outer radius to the 
inner radius, $ =r, - r;,— wall thickness. 

The circumferential thermal stress o,,,, on the inner 
surface is calculated as follows 

2 
ee a (8) 
l-v 

The circumferential thermal stress on the outer surface 

is calculated from the similar formula 


2 
= FE MS (9) 
l-v K 
where the so-called shape factors for the internal sur- 
face ¢,,, and the external surface ¢,, are given by the 
formulas 


(7) 


Oo. 


ww 


WZ 


(10) 


(u? -1) (3u? -1)-4u* Inu 
(2 —— ass 


(u ~1)(u-1y’ 
1 (4 ~1)-4u? Inu 


8 (u? -1)(u -1) 

Similarly, the average temperature over the thickness 
of the sphere can be derived from Eq. (6). The shape 
factors for a sphere: ¢,, for the inner surface and Ø, for 
the external surface, are given by 


1 9u° —5u° — 5u? +1 

Prw — 15 G -1)(u iy (12) 

T 1 (už -1)-5u? (u-1) 
10 (u° —1)(u-1)" 


The shape factors at the inner and outer surface of cy- 
lindrical and spherical walls are presented in Fig. 3. 

The absolute values of the shape factors for the 
spherical wall øw and øx are smaller than the appropriate 
values of the shape factors for the cylindrical wall ¢,,,, or 


uz (Fig. 3). 
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Fig. 3 Shape factors for the inner and outer surface of cylin- 
drical and spherical walls. 


For the same rate of temperature changes, the thermal 
stresses in the spherical walls will be greater than the 
thermal stresses in the cylindrical walls. 


Optimum Heating of Components Without Openings 
Heating or cooling of the components is conducted in 

a way, which will assure that the circumferential thermal 

stress on the inner surface equals the allowable stress. 


Oww =O, (14) 
Substituting Eq. (8) into Eq. (14) gives 
2 
BES 4, = 0, (15) 
l-v 


After transformation of Eq. (15), the expressions for 
the heating or cooling rate of the thick walled, cylindrical 
or spherical boiler components are obtained 

cylindrical wall 

O,K 


T (16) 
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The temperature change rate is positive during the 
heating process and negative during the cooling process. 
The wall temperature depends on the fluid temperature, 
but the temperature difference between the selected 
points in the components is constant and independent of 
time and the value of the fluid temperature. 


Quasi — Steady Distributions of the Temperature and 
Stresses in Pressure Components with Openings 

The quasi - steady state temperature distribution in 
components of complex shape can be determined using 
Eq. (1) with the appropriate boundary conditions. Equa- 
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tion (1) describes a steady state temperature distribution 
in the body T(zt) with uniformly distributed heat sinks 
with the specific power equal to: ġ, = v, / K. 

In the following, an application of the quasi — steady 
state approach for determining allowable heating rates in 
components with complex shape will be discussed. Since 
stresses at points P; and P, (Fig. 4) limit the heating rate 
of the boiler drum or other thick-walled vessels weak- 
ened by nozzles [3-6], thermal and pressure stresses are 
calculated for these points. 


Fig. 4 Section of the boiler drum — downcomer junction and 
mesh of finite elements used for thermo-mechanical 
analysis of the drum — downcomer junction 


The following geometrical dimensions of the boiler 
drum: inner diameter of the drum D = 1.7m, inner di- 
ameter of the downcomer d = 0.09m, thickness of the 
drum wall H = 0.09 m, thickness of the downcomer wall 


h = 0.006 m and physical properties of the boiler material: 


k= 42 Wi(mK), c = 538.5 I/(kg-K), p = 7800 kg/m’, E = 
1.96x10!! N/m’, 8 = 1.32x10° 1/K and v= 0.3 were as- 
sumed for the calculations. At first, time changes of the 
circumferential stress in time at point P, for different heat 
transfer coefficients were computed (Fig. 5). An analysis 
of the results shows that the quasi — steady state at low 
values of the heat transfer coefficient h = 500W/(m’-K) is 
set only after a long heating time at a constant rate. After 
4000 seconds of heating the stresses did not reach the 
constant values, characteristic for the quasi — steady state. 
Equivalent stress distributions in the boiler drum — 
downcomer junction determined using ANSYS for dif- 
ferent pressures and heating rates are depicted in Figures 
6, 7, and 8. The largest stress due to the pressure is at the 
point P; which is located at the hole edge (Fig. 6). 

The equivalent stress at the point P; is about five times 
larger than the stress at the point P) when the drum is 
loaded only by the inner pressure. If the pressure inside 
the drum is equal to zero, then the equivalent thermal 
stress at the point P, is higher in comparison with the 
stress at the point P; (Figs 7a and 7b). When pressure is 
increased to design pressure p = 10.87 MPa then the 
stresses due to pressure for vr = 1 K/min at the point P, 
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exceed thermal stresses and the largest equivalent stress 
occurs at the point P, (Fig.8a). After increasing the 
heating rate to vr = 12 K/min negative thermal stresses 
exceed positive stresses due to pressure and the largest 
equivalent stress is at point P» (Fig. 8b). 


v = 3 K/min H 
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Fig. 5 Time variations of the circumferential stress o, in time 
at point P, for different heat transfer coefficients A. 
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Fig. 6 Equivalent stress distribution caused by inner pressure 


The stress concentration coefficients a at the point P, 
and P, (Fig. 4), for the quasi — steady state as a function 
of the heat transfer coefficient are presented in Fig. 9. 
Thermal stresses determined from Eqs (8) and (9) refer to 
simple shapes, such as hollow cylinders and spheres. In 
real life, components used in power units often have 
openings, which can be viewed as material discontinuity 
causing stress concentration, what shall be taken into 
consideration when determining the allowable heating or 
cooling rate. Most often the assumed value of the thermal 
stress concentration coefficient at the edge of the opening 
is: @r= 2. Then, the influence of the heat transfer coeffi- 
cient A on the thermal stress concentration coefficient ær 
at the edge of the opening is not considered. 
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a) 


0.018206 8.79698 17.5757 26.3545 35.1333 
x ~—S Ee a 


4.40759 13.1864 21.9651 30.7439 39.5227 


b) 
21.488 79.0904 136.693 194.295 251.897 


50.2892 107.892 165.494 223.096 280.699 


Fig. 7 Equivalent stress distributions at the drum — down- 
comer junction without pressure load (gauge pressure 
p =0 MPa) for different heating rates; a) vr = 1 K/min, 
b) vr = 12 K/min 


a) 


0.218469 105.564 210.909 316.254 421.599 


52.8911 158.236 263.582 368.927 474.272 


b) 


5.74966 96.0498 186.35 276.65 366.95 
eS 
50.8997 141.2 231.5 321.8 412.1 


Fig. 8 Equivalent stress distributions at the drum — down- 
comer junction without pressure load (p = 10.87 MPa) 
for different heating rates; a) vr = 1 K/min, b) vy = 12 
K/min 
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Fig. 9 Thermal stress concentration coefficient ær as a func- 
tion of the heat transfer coefficient h 


However, the thermal stress concentration coefficient 
at the opening edge in the complex shaped components 
depends on the heat transfer coefficient h. Since the 
boiler drum has openings for downcomers, an intense 
concentration of stress is present at the edge of the open- 
ing. It can also be noted that the stress concentration co- 
efficient at the edge of the opening strongly influences 
the heat transfer coefficient on the inner surface of the 
boiler drum and the downcomer. The value of the stress 
concentration coefficient a; is lower at point P; — as 
compared to its value at point P, — regardless of the value 
of the heat transfer coefficient (Fig. 9). 

When the heat transfer coefficient is equal to: h = 
10000 W/m’K, the stress concentration coefficient at 
point P; is: ær = 1.53, and at point P3: ær = 1.86. The 
thermal stress concentration coefficient was calculated 
using the formula 


ap =—— (18) 
Cis, 


where o, p denotes circumferential stress at the edge 


of the opening at point P|, and o,,,=0,(7,) stands 


o 
for circumferential stress on the inner surface of the boi- 
ler drum or pressure vessel without openings. 

The total circumferential and equivalent stress at the 
edge of the opening at points P; and P, can be approxi- 
mated linearly as a function of the heating rate. At the 
beginning of the boiler start — up, the gauge pressure in 
the drum is equal to zero thus the absolute values of 
circumferential and equivalent stresses at the points Pı 
and P, are equal (Fig. 10). The total circumferential and 
equivalent stresses presented in Fig.11 were computed 
considering pressure and thermal loads. 

The circumferential stress at points P; and P, for pres- 


sure p = 10.87 MPa can be approximated by a linear 
function with respect to vr 


O, =a+bvr (19) 
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Fig. 10 Circumferential stress and equivalent stress at the 
edge of the downcomer opening at points P; and P, 
for the gauge pressure p = 0 MPa; 1 — circumfere- 
ntial stress at point P;, 2 — equivalent stress at point 
P,, 3 - circumferential stress at point P), 4 — equiva- 
lent stress at point P», 5 — allowable circumferential 
stress, 6 — allowable equivalent stress 
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Fig. 11 Circumferential stress and equivalent stress at the 
edge of the downcomer opening at points Pı and P3 
for the pressure p = 10.87 MPa; 1 — circumferential 
stress at point P4, 2 — equivalent stress at point P4, 3 
circumferential stress at point Po, 4 — equivalent 
stress at point P), 5 — allowable circumferential stress, 


6 — allowable equivalent stress 


The allowable vr; and vr: boiler drum heating rates 
can be determined from the condition o= oa, where the 
circumferential stress o, is given by Eq. (19). The same 
method can be applied to other components of complex 


shapes. The issue of optimization of heating and cooling 
of thick walled boiler components was presented in [5] 
without the assumption of the quasi-steady state. 

However, the developed method for determining the 
allowable rate of cooling or heating pressure components 
is simple and can be easily used in the practice. 


Components with Non - uniform Circumferential 
Temperature Distribution 


There is an uneven temperature distribution on the 
drum circumference due to the presence of the water and 
steam zone. During start-up of the boiler, steam con- 
denses on the inner surface of the drum in the steam zone. 
During the steam condensation, the upper part of the 
drum is heated very quickly because the heat transfer 
coefficient at the time of condensation is very high. The 
bottom of the drum, occupied by the water, the coeffi- 
cient of heat transfer from the water to the inner surface 
of the drum is low, and the heating of the lower part of 
the drum occurs much more slowly. Heating at different 
rates of the upper and lower portions of the drum results 
in a temperature difference between the top and bottom 
of the drum. It should be pointed out that the temperature 
of the upper part of the drum is higher than the tempera- 
ture of the lower part. A similar phenomenon can be 
observed in the horizontal sections of the pipeline and 
superheater headers. At the beginning of the start-up of 
the boiler, pipes are cold and water vapor condenses on 
the walls of the pipes. At the bottom of the horizontal 
piping and other pressure components of the boiler water 
accumulates at the bottom, especially when the drainages 
of the elements are not functioning properly. Such a case 
occurs when the horizontal pipelines and superheater 
headers are plastically deformed. In a later stage of the 
boiler start-up superheated steam flows through the upper 
part of the pipeline while at the bottom of the pipeline is 
saturated water. This is the reason for the occurrence of 
very large temperature differences on the circumference 
of the pipe, which reach up to about 200K. For this rea- 
son, permanent deformation of pipes and headers often 
occur, which in turn makes that the draining horizontal 
elements are not working properly. The stresses caused 
by temperature non-uniformity on the circumference of 
the pipe may be higher than the stress caused by temper- 
ature difference across the thickness of the pipeline. In 
the current European Standard EN 12952-3, thermal 
stresses caused by the temperature difference on the cir- 
cumference of the cylindrical element are not at all taken 
into account. In this paper, an analysis of the heating of 
horizontal thick walled components will be carried out 
taking into account the different values of heat transfer 
coefficients in the water and steam zone (Fig. 12). Cal- 
culations are performed assuming that the temperature of 


the fluid increases linearly with time so that after some 
time forms a quasi-steady state temperature field. This 
paper examines thermal and strength loads of the drum of 
the steam boiler OP-210M during startup. Figure 13 illu- 
strates the time changes of the water level. The water 
level in the drum situated H„ = 10 cm below the hori- 
zontal plane passing through the axis of the drum is taken 
as the zero level. The water temperature in the evaporator 
increases resulting in an increase in the specific volume 
of the water and rising water level in the drum. After 
about 12000 seconds, the excess water is discharged from 
the drum and the water level is lowered to a level close to 
zero. The time variations of the saturation pressure in the 
drum is depicted in Fig. 14. Measured time variations of 
the drum temperature at seven points shown in Fig. 14 
indicate that temperature differences on the drum cir- 
cumference, especially at the initial period of the drum 
heating are significant. 


The cross-section of the boiler drum with locations of 
thermocouples used for temperature measurement 


Fig. 12 
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Time variations of the level of water in the drum of 
OP-210M boiler during the start-up 


Fig. 13 


As already explained earlier, the faster heating of the 
upper part of the drum is due to vapor condensation on 
the inner surface of the drum, during which the heat 
transfer coefficient is very high. After about 15000 
seconds the temperature difference between the top and 


the bottom of the drum is reduced considerably (Fig. 14). 
The maximum temperature difference between the top 
and bottom of the drum does not exceed 50 K. It must be 
emphasized that such a large temperature difference at 
the circumference of the drum is the cause of high ther- 
mal stresses. The analysis of the experimental results 
shows that in the initial stage of drum heating there are 
large temperature differences on the perimeter, which are 
the cause of additional thermal stresses. The current boi- 
ler standards do not take into account the boiler stresses 
due to temperature difference on the perimeter of the 
horizontal elements. Only quasi-stationary thermal 
stresses caused by the temperature difference across the 
wall thickness are considered. 

Assuming the existence of quasi-stationary tempera- 
ture field in the pressure element, the temperature differ- 
ence on the component circumference caused by different 
heat transfer coefficients on the inner surface will be 
taken into account. The time variations of the circumfe- 
rential stress at the edge of the drum — downcomer junc- 
tion at the points P; and P, are shown in Figs. 15 and 16. 
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Fig. 14 Temperature and pressure variations on the drum 
circumference during the boiler start —up 


The analysis of the results depicted in Figs. 15 and 16 
shows that the absolute value of thermal stresses at the 
point P, is much greater than those at the point Pj. It 
should be noted that the heat transfer coefficient in the 
water area is greater than 1000 W/(m°-K). Due to the 
higher value of compressive thermal stress and a lower 
value of the tensile stress due to pressure occurring at 
point P5, the permissible rate of the drum heating should 
be determined due to the stress at the point P2. Changes 
in the circumferential stress and equivalent at points P, 
and P, as a function of heating rate are shown in Figure 
17. In Figure 17a, the gauge pressure in the drum is p = 0, 
while in Figure 17b the pressure is equal to the design 
value of 10.87 MPa. 

It should be emphasized that during the start-up of the 
boiler heat transfer coefficient in the water zone is less 
than in the steam zone. In such cases, a maximum al- 
lowable heating rate is determined by the stresses at the 
point P}. 
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Fig. 15 Circumferential stresses on the hole edge at the point 
Py 
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Fig. 16 Circumferential stresses on the hole edge at the point 
Py 


For a typical case that occurs in practice [7], when h,, 
= 1000 W/(m°K) and h, = 10000 W/(m*K), the permiss- 
ible heating rate at the beginning of the boiler start-up 
when p = 0 is approximately 3 K / min. Currently, the 
permissible rate of the drum heating recommended by 
boiler manufacturers is about 1 K/min. Allowable heating 
rate equal 3K/min (Figure 17a) resulting from the FEM 
analysis is significantly higher than the value prescribed 
by the boiler manufacturers. The start-up of the boiler 
can thus be shortened three times. 


Conclusions 


Summarizing the obtained results, it should be noted 
that the assumption of quasi — steady temperature field 
significantly simplifies the determination of the allow- 
able rates of fluid temperature changes during the heating 
and cooling of boiler thick walled components. The 
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thermal stress concentration coefficients can be quickly 
determined with the use of the Finite Element Method 
(FEM), also for components of very complex shapes. For 
flat, cylindrical and spherical walls, in which the tem- 
perature distribution is one-dimensional, the thermal 
stress concentration coefficient does not depend on the 
value of the heat transfer coefficient. It was also shown 
that the heat stress concentration coefficient for the 
structural components of complex shapes depends on the 
value of the heat transfer coefficient on their inner sur- 
face. 
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Fig. 17 Circumferential stress and equivalent stress at the 
edge of the downcomer opening at points P; and P, 
for hy, = 1000 W/(m"K) and h, = 10000 W/(m*K); a) 
drum gauge pressure p = 0 MPa, b) drum gauge 
pressure p = 10.87 MPa; | — circumferential stress at 
point Pı, 2 — equivalent stress at point P}, 3 — 
circumferential stress at point P), 4 — equivalent 
stress at point P), 5 — allowable circumferential 
stress, 6 — allowable equivalent stress 


It has also analyzed quasi-steady state temperature 
field and stress in the drum — downcomer junction con- 


sidering the different heat transfer coefficients in the wa- 
ter and steam space. Graphs presented in the chapter al- 
low to determine heating rate limits at the beginning of 
the boiler start-up when the pressure is equal to 0 and at 
the end of starting when the pressure is equal to the 
nominal pressure. Permissible drum heating rates, as- 
suming both circularly symmetric and asymmetric circu- 
lar cross-section temperature field in the due drum ob- 
tained by the proposed method are several times higher 
than the limit values recommended by the manufacturers 
of boilers. 
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